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Abstract. We shortly recall the derivation of the Faddeev-Yakubovsky differential equa- 
tions and point out their main advantages. Then we give a review of the numerical ap- 
proaches used to solve the bound- state and scattering problems for the three- and four-body 
systems based on these equations. A particular attention is payed to the latest developments. 



1. Brief history 

A crucial step in the theory of few-body quantum systems has been done in 1960 when 
Faddeev introduced his celebrated integral equations O. The Faddeev equations became 
a basis for the three-body scattering theory with short-range interactions. They have also 
been used through the years in numerical calculations. As for the Af-body systems with 
N >3, Faddeev's idea of an explicit cluster-channel separation was completely elaborated by 
1967 when Yakubovsky's integral equations became available Q. The Faddeev differential 
equations, in their S- wave two-dimensional version, were first discovered in 1968 by Noyes 
and Fiedeldey Q. The Faddeev differential equations in their complete form have been 
analysed in 1976 by Merkuriev, Gignoux and Laverne [4j. In particular, they found and 
approved the asymptotic boundary conditions that are needed to be added to the equations in 
order to find unique physical solutions corresponding to various scattering processes. 

It was the great success of the Faddeev differential equations that stimulated Merkuriev 
to look for a generalization of these equations to few-body systems with arbitrary number of 
particles. By 1982 he has successfully solved the problem jointly with his then Ph.D. student 
Sergei Yakovlev 0. They not only found the "right" differential equations for the Yaku- 
bovsky components but also described the boundary- value problems for these equations that 
correspond to certain scattering processes L6j. In a somewhat more abstract form, the Ya- 
kubovsky differential equations can be seen in paper [7j by Benoist-Gueutal and L'Huillier, 
published also in 1982. 

In Section 2 we recall how the Faddeev-Yakubovsky differential equations look like. The 
best way to do this is simply to derive them. We concentrate only on the algebraic context of 
the derivation, leaving apart the description of the scattering boundary conditions. In Section 
3 we review the numerical approaches used to solve the Faddeev-Yakubovsky differential 
equations. A particular attention is payed to the latest developments. Finally, in Section 4 
we mention some still challenging problems. 
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2. FADDEEV AND YAKUBOVSKY DIFFERENTIAL EQUATIONS 

Let Ho be a linear (not necessarily Hermitian) operator on a Hilbert space $), and V another 
linear operator on $). At the moment we specify neither Ho, nor V. These operators may be 
of completely arbitrary nature. Temporarily assuming that both Ho and V are bounded we 
avoid complications with domains of the operators involved and concentrate mainly on the 
algebraic part of the Faddeev-Yakubovsky scheme. The only essential assumption is that by 

some reason the perturbation V is split into the sum V = V\ + V% H V V n of n (2 < n < oo) 

terms V a . Surely, one is interested in the perturbed operator H = Ho + V. 

Suppose that z is an eigenvalue of H and *P the corresponding eigenvector. Assume, in 
addition, that z does not belong to the spectrum o(Ho) of the unperturbed operator Ho. Then 
Ho — z is invertible and (Ho + V)^ = z^ is equivalent to the Lippmann-Schwinger equation 

y = -(H -z)- 1 £ V a y. (1) 
a=\ 

Introduce the vectors 

¥a = -(Ho ~ Z)~ l VaW, a = 1,2, . . . ,n. (2) 
These vectors are called the Faddeev components of the eigenvector X V. Obviously, combin- 

n 

ing (OD) and © implies Wb = ^« Taking this into account, rewrite © in the form 

J3=i 

Then apply to both sides of © the operator Ho~z and arrive at the equations 

(#o + V a -z)y/ a = -V a £ Vj3 . (4) 

These are just the celebrated Faddeev "differential" equations for the components yr a of 
an eigenvector of H. A remark on the spectrum of the associated Faddeev operator can be 
found in Appendix. 

For z o(Ho + Va), oc = 1,2, . . . ,n, one can invert the operators Ho + V a — z and, thus, 
rewrite © in the equivalent form 

^ = -(^0 + V a -z)- 1 Va £ Yp. (5) 

Equations © are known as the Faddeev integral equations. 

The Faddeev equations © turn into truly differential equations if //o is a differential op- 
erator as this happens in the case of few-body problems with pairwise interactions, studied 
in coordinate representation. In this case Ho is simply the kinetic energy operator in the 
center-of-mass frame and V a 's stand for the two-body potentials. 

In the three-body case the Faddeev equations © or © represent simultaneously the first 
and the last step in the Faddeev-Yakubovsky approach. This reflects the fact that with suffi- 
ciently rapidly decreasing and smooth potentials V a , oc = 1, 2, 3, the four times iterated Fad- 
deev integral equations © are Fredholm and even compact (unlike the Lippman-Schwinger 
equation © that remains non-Fredholm after any number of iterations, even for the complex 
energies z). 
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Now we turn to the four-body case. It is convenient to depict the four-body system by 
an icon like IxJ where the dots denote the particles and the connecting line segments are 
associated with both the particle pairs and respective two-body interactions. There are n = 6 
different pairs (and potentials V a ) that are numbered by a = 1,2, . . . , 6. 

It turns out that the four-body Faddeev integral equations ® are not Fredholm even 
for complex z. This can be understood already from the fact that in the four-body coordinate 
space R 9 the supports of any two finite non-zero two-body potentials V a and Vp, j8 ^ a, have 
an unbounded cylindrical intersection. 

Yakubovsky's recipe to tackle the four-body problem consists in the following. First, 
introduce the two-cluster partitions: ]/! ?\I 'A I I IT X . Number them by the Ro- 
man letter a (or b, c etc.). Then say that a pair a belongs to a partition a and write a C a 
if the two-body subsystem a belongs to one of the clusters in the partition a. For example, 
I o pi 11 while \ c 57 t£ . In such a case the sequence aa is called the chain 

of (consecutive) partitions. It is assumed that the partition a consists of one two-particle 
cluster formed by the pair a and two one-particle clusters obtained from the partition a by 
breaking the corresponding "bonds". It is easy to see that there are 18 different chains of 
partitions: 12 due to 3 + 1 and 6 due 2 + 2 starting partitions. 

Using the four-body Faddeev integral equations © introduce the following new vectors, 
called the Yakubovsky components (of the eigenvector *P): 



Vaa = -(Ho + V a -z) l V a £ V/3- (6) 



Then apply to both parts of © the operator Hq + V a — z and obtain 



(H + V a -Z)\lfaa = ~V a £ Vp = -V a £ £ V*j8 ( 7 ) 



taking into account that Y,bDfi Vfib = Yp- Further, by applying some combinatorics one ar- 
rives at the identity £ £ \j/ b p = £ £ \j/ b p and then from © it follows that 

(P^a)ca bDfi b (Pjta)ca 



(H + V a -z)\l/aa+V a £ Wafi = ~V a £ £ ( 8 ) 



These 18 equations for the 18 unknowns are just the desired Yakubovsky differential 
equations. 
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The table above demonstrates the remarkable structure of the 18 x 18 block operator ma- 
trix associated with the Yakubovsky equations ([8]). The icons of three-cluster partitions depict 
the indices a of the potentials V a in the non-zero entries of the matrix. 

3. Numerical approaches, ideas and tricks 

Compared to the Schrodinger equation, the main numerical advantage of the Faddeev- 
Yakubovsky differential equations is that the asymptotical boundary conditions for their 
physical solutions imposed at infinitely large distances are much simpler than those for the 
total wave function, and they can be much easier incorporated into the numerical scheme, 
especially in the case of scattering processes. Compared to the Faddeev- Yakubovsky inte- 
gral equations, it is much easier to tackle the local two-body potentials. Typically, with such 
potentials the matrices of the discretized differential equations have a band structure unlike 
in the case of the integral ones. 

If particles are identical, the number of the essentially different Faddeev- Yakubovsky 
components reduces. The components corresponding to different indices may be obtained 
from each other by a simple unitary transformation. For example, the Yakubovsky equations 
([8]) for four identical bosons reduce to a system of only two equations (8|. 

The next step is to reduce the dimension of the equations by making a partial-wave de- 
composition. Since the total angular momentum L is a conserving quantity, projections onto 
the reducing subspaces associated with a fixed value of L keep the emerging equations exact. 
At this stage the three Eulerian angles fixing the particles plane turn out to be separated and 
the dimension of the Faddeev equations reduces to 3. This is a modern development to stop 
at this stage and solve the three-dimensional Faddeev equations. Such an approach was first 
developed in 19 ]. It is especially well suited to the case where two-body interactions depend 
only on the distance between particles. The three-dimensional Faddeev differential equations 
were first solved numerically in fLOl . 
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But historically, the first were the two-dimensional partial-wave (integro-) differential 
Faddeev equations where, for a fixed value of L, another variable, the angle between Ja- 
cobian vectors, was eliminated ifTTTl . Further reduction, already to a one-dimensional form, 
consists in performing the hyperspherical adiabatic expansions (see, e.g., fL2lO . 

As for the discretization of the equations with respect to the spatial variables, the finite- 
difference approximation was initially used [llj. Starting from Q31, for this purpose various 
spline approximations are often employed. An interesting trick with a tensor factorization 
of the left-hand sides of the discretized Faddeev equations has been suggested in fL4ll . The 
factorization trick works very well within iterative approaches (see fl31l . |[T6lO . 

Four-body calculations based on the Yakubovsky differential equations are still rather 
scarce. Recent examples of such calculations can be found in |[T7l[T8l . 

Latest developments in approaches to solving the Faddeev- Yakubovsky differential equa- 
tions are related to the case where the interaction potentials have an extremely strong re- 
pulsive component at small distances between particles (like, e.g., in atom-atom potentials). 
One of the options to tackle such potentials is in approximating the repulsive part by a hard 
core. A generalization of the three-body Faddeev differential equations to the hard-core in- 
teractions has been done yet in 1983 by Merkuriev and the present author lfT9ll . Differential 
Yakubovsky equations for this model were derived in [20 j. According to HUE!, the effect 
of hard cores is reproduced by imposing two-sided boundary conditions on the Faddeev- 
Yakubovsky components at x a = c a , a = 1 , 2, . . . , where x a denotes the distance between 
particles of the pair a and c a the sum of their core radii. In the three-body case these condi- 
tions are quite natural: Y?p=\ Wp \ x =c = 0, a = 1,2,3, while in the four-body one they are 
more nontrivial: 



It is assumed that V a = for x a < c a and that equations © and ([8]) are considered on the 
respective whole coordinate space R 6 or R 9 , except the hypercylinders x a = c a where the 
hard-core boundary conditions are imposed. Thus, when a potential with a strong repulsive 
core is replaced by the hard-core model, one approximates inside the core domains only the 
kinetic energy operator Hq instead of the sum of Hq and the huge repulsive term. In this way 
a much better numerical approximation is achieved. 

In recent years this approach was extensively used to calculate binding energies, reso- 
nances, and scattering observables for the systems of three helium atoms (see papers ll2Tll22ll 
and references therein). In 2006 the Yakubovsky differential equations ([8]) with the hard-core 
boundary conditions © were employed for the first time and with a big success by Lazaus- 
kas and Carbonell lf23ll to calculate properties of the four-atomic 4 He4 system. They have 
also performed 4 He3 calculations with the hard-core Faddeev differential equations. 

Another (purely numerical) approach to tackle the two-body potentials with extremely 
strong repulsion at small distances has been developed by Roudnev and Yakovlev. A quite 
detail description of this approach is given in lfl5ll . 



For years the Faddeev- Yakubovsky differential equations showed their high suitability 
and efficiency when one solves the three- and four-body problems, particularly the scattering 
ones. Recent advances are related to solving the problems with hard-core (or practically 




for any chain aa. 



(9) 



4. Challenges 
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hard-core) interactions. Still challenging, however, is the numerical solving of the scattering 
problems that involve triple collisions. Also many questions remain unanswered in the case 
of few-body scattering problems with Coulomb interactions, especially with the attractive 
ones. 

Appendix: On the spectrum of the Faddeev operator 
Faddeev equations © represent the spectral problem for the n x n block operator matrix 

/flb + Vi Vi ... Vi \ 

„ V 2 H + V 2 ... V 2 

H F = 

V V n V n ... Ho + V n J 

n times 

~ , A s 

considered as an operator on the Hilbert space S) = # 5} ... 5}. The matrix Hp is called the 

Faddeev operator associated with a Hamiltonian of the form H = Ho + V\ + V 2 H h V n . 

The passage from H to Hp adds to some auxiliary spectrum. But this auxiliary spectrum is 

nothing but that of Ho. More precisely, o(Hf) = o(H) U o(Ho) (for a proof see, e.g., l24l ). The spec- 
trum o(Ho) is usually called the spurious spectrum of the Faddeev operator, although it is quite not 
"dangerous", since it is assumed that one knows everything about the unperturbed Hamiltonian Ho. 
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